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Hydrogenation of graphene on the a-quartz (0001) SiCb substrate is studied, considering different 
surface terminations in order to take into account the amorphic nature of the material. Our ab 
initio calculations show that the formation of graphane by hydrogen adsorption on graphene is 
energetically favored on hydroxy 1 and oxygen terminated surfaces, whereas silicon termination and 
reconstruction of the oxygen termination hinder adsorption. Our results indicate that in order to 
fabricate graphane on Si02, it is beneficial to oxygenize the surface and saturate it with hydrogen. 
For the pristine graphene on the substrate, we find only marginal changes in the low-energy band 
structure for all surface terminations. 



I. INTRODUCTION 

The two-dimensional allotrope of carbon, graphene, 
has many unique properties. For example, the disper- 
sion is linear at the Fermi level, which makes the mate- 
rial semimetallic. For some applications like field-effect 
transistors, a band gap would be more desirable. Luck- 
ily, there are ways to induce a gap, like cutting graphene 
into narrow nanoribbons. 1 An alternative approach is the 
chemical functionalization. In particular, hydrogen ad- 
sorption is an interesting option for two reasons: firstly, 
a band gap is induced 2 -* 3 - and secondly, the hydrogenated 
graphene sheet could serve as a hydrogen storage^ 

In experiments, the graphene layer is typically placed 
on a substrate, such as silicon dioxide, 5,6 silicon 
carbide^ - — or a metal surface. 2,10 The layer is then ex- 
posed to hydrogen by using hydrogen plasma 5 or atomic 
hydrogen^ - — The dense structure of the carbon layer 
suggests that hydrogenation should occur only on one 
side of the carbon planed The hydrogen coverage reached 
in experiments varies and is difficult to estimate. Sessi et 
all demonstrated a quite dense and even hydrogen ad- 
sorption on SiC. Guisinger et aim reached only a "very 
low" hydrogen coverage and Balog et al. — were not able 
to estimate the coverage on the same substrate. On Si02, 
neither Elias et all nor Ryu et al. 6 report the coverage. 
On the metal substrates, Haberer et alll find a maxi- 
mum hydrogen coverage of 25% and Balog et al? report 
3-27%. The band gap of substrate-deposited graphane 
has not been accurately measured but a lower limit of ap- 
proximately 0.5 eV has been estimated 2,5 . Interestingly, 
in addition to hydrogenation, Sessi et al. 7 also demon- 
strated selective dehydrogenation of the structure that 
restores the properties of pristine graphene. 

This far, theoretical studies on hydrogenated graphene, 
have concentrated on the suspended layer 
Graphane, with a hydrogen atom attached to every car- 
bon atom, is found to be a good insulator with a large en- 
ergy gap of even 6 eV in theoretical calculations ! 20 ! 21 In- 
stead, a freestanding graphene membrane hydrogenated 
only on one side (called also graphone 12 ), has not been 
found stable. 12,13 This opens the question of the role of 
the substrate for the structural and electronic properties 



of graphane. 

In theoretical calculations, the presence of a sub- 
strate is widely neglected even in the case of pristine 
graphene. Additionally, the theoretical descriptions^ 2 -"— 
of graphene on the Si02 substrate are contradictory even 
as to whether or not graphene forms bonds to the sub- 
strate, destroying the linear dispersion at the corners of 
the Brillouin zone. To the best of our knowledge, the 
hydrogenation of graphene on Si02 has not yet been ad- 
dressed. 

In this work, we study the hydrogenation of graphene 
on Si02 using ab initio calculations. We consider several 
different surface terminations in order to yield insight 
into the hydrogenation process of the real amorphic Si02- 
Besides the simply cleaved oxygen terminated (OT) and 
cleaved Si terminated (SiT) surfaces, we also consider 
the reconstructed oxygen terminated (ROT) and hydro- 
gen saturated oxygen terminated (OHT) surfaces, since 
the OT surface is reactive and prone to reconstruction or 
hydrogen saturation in ambient conditions^"— We con- 
sider hydrogen coverages of graphene ranging from 1/8 
to 3/4 monolayers (ML). Our calculations show that the 
hydrogen atoms are most likely to attach onto graphene 
on the OHT and OT surfaces. On the OHT surface, the 
arrangement of the hydrogen atoms in the lowest-energy 
configuration resembles the chair configuration of free- 
standing graphane. On the OT surface, hydrogen atoms 
form lines and make structure dispersive mainly in one 
direction. On the ROT or SiT substrates, we do not find 
stable hydrogen configurations. Our results suggest that 
in order to facilitate hydrogen adsorption on graphene, 
the substrate should be oxygenated and hydrogen passi- 
vated. 



II. METHODS 

Our calculations were performed using density- 
functional theory (DFT) with a van der Waals (vdW) 
correction, 30 implemented in the all-electron FHI-aims 
code^i Double numeric plus polarization basis set of nu- 
merical atom-centered orbital basis functions and PBE 32 
exchange correlation functionals were used. The struc- 
tural relaxation was converged until the forces acting on 



2 



the atoms were less than 0.001 eV/A and the electronic 
degrees of freedom were convert ed below 10~ 6 eV. Due 
to the abundance of met ast able intermediate states, tight 
convergence criterions were found to be important. In the 
relaxations, no spin polarization was taken into account. 
For some of the final geometries, however, we searched 
spin moments but did not find any spin polarization. Of 
course, this does not preclude the possibility of more com- 
plicated geometries with finite spin moments. 

Silicon dioxide in the a phase was modelled using three 
Si02 unit cell thick slabs, corresponding to a layer of 
15.8 A. Periodic boundaries were applied and a vacuum 
layer of approximately 20 A was placed between the pe- 
riodic images of the slabs. Four different Si02 surfaces 
were created: SiT, OT, ROT, and OHT surfaces. The 
starting coordinates for the relaxations for the ROT sur- 
face were taken from the work of Goumans et al*2! The 
lowest-energy OHT surface was obtained by relaxing the 
OT surface in the presence of additional hydrogen atoms, 
using different initial positions for the surface hydrogens. 

The graphene layer was then placed on both sides of 
the silicon dioxide slab, the two interfaces being equiva- 
lent but rotated by 60°. The lattice mismatch between 
the unit vectors of graphene and Si02 was only 1.3 %. In 
each Si02 unit cell, there were four graphene unit cells 
corresponding to eight carbon atoms (2x2 carbon super- 
cell). A 6x6 /c-point grid was used for this supercell. 
The carbon atoms and three uppermost substrate atom 
layers were relaxed. Different initial positions for the car- 
bon grid were used and they were all found to relax to 
the same state. 



III. GRAPHENE ON Si0 2 SURFACE 

For the relaxed geometries, we find only weak bind- 
ing of graphene on all surfaces. Throughout this paper, 
with strong binding we refer to structures with substrate- 
carbon bonds, as opposed to weak binding characterized 
by longer substrate-carbon distances. The graphene layer 
is farther away from the substrate for the SiT surface 
(3.4 A), compared to the other surfaces («3 A). 33 The 
relative position of the graphene layer on the four differ- 
ent surfaces is presented in Fig.[TJ For the SiT, ROT, and 
OHT surfaces, we get the hollow adsorption side, which 
fully agrees with previous work* 2 ^ - — For the OT surface, 
the uppermost oxygen atom of the substrate is slightly 
displaced from a bridge position between two carbons 
towards the center of the hexagon. Previous simulations 
have suggested top 2 ^ 2 ^ or bridge 26 positions. 

Otherwise, previous DFT studies have not found agree- 
ment on whether graphene strongly binds to the Si02 
surface ^ 22 ! 24 or not j 23 i 26 This discrepancy can also be seen 
in the values for the equilibrium distance of graphene 
from the substrate that range from 1.43 A 22 to approx- 
imately 3.4 A. 24 During our relaxations, we find inter- 
mediate geometries with small forces that still continue 
to relax to clearly lower-energy geometries through slow 



sliding of the carbon layer on the surface. Because some 
of the intermediate states were bound to the surface, one 
possible cause for the discrepancy is the effect of numer- 
ical settings in the different calculations. At the very 
least, the relaxations need to be done with a high nu- 
merical accuracy, relaxing the forces until they are very 
small. Also the inclusion of the van der Waals correction 
can, in principle, affect the results but it should provide 
only additional attraction in comparison to the non-vdW 
case. 

The opening of the gap is connected to how tightly 
the graphene is bound to surface. In our calculations 
the gaps are zero, apart from the OT surface where 
we find a gap of approximately 10 meV. The band 
structure of graphene on all different surfaces shows 
the linear dispersion near the K-point, characteristic 
for freestanding graphene (see Fig. [TJ). Some substrate 
bands cross the graphene bands but only on the OT 
and OHT surfaces they are near the Fermi energy. In 
the previous calculations on the OT surface, either 
the linear dispersion without a gap is preserve d 23 ! 26 or 
the graphene layer binds strongly, leading to a large 
gap j 22 i 24 On the ROT and OHT surfaces, a small gap 
of approximately 20 meV has been reported* 2 ^ In the 
experiments, the linear dispersion has been observed for 
Si02-deposited graphene^ 



IV. GRAPHANE ON Si0 2 SURFACE 

When hydrogen atoms adsorb to the carbon sheet, the 
conducting graphene is transformed into a semiconduct- 
ing graphane. In the freestanding fully hydrogenated 
graphane sheet, hydrogens attach to the carbons of dif- 
ferent sublattices on different sides of the sheet. The 
estimates for the semiconducting energy gap depend on 
the level of theory, ranging from approximately 3 eV of 
DFT to 6 eV in GW calculations ^OLi&so The carbon- 
carbon bond length in graphane, 1.5-1.6 A ^ 12 i 15 i 18 ~ — is 
longer than in graphene, 1.42 A. The longer bonds are 
counterbalanced by the non-planar structure of the car- 
bon layer. Using a pristine graphene sheet and molecular 
hydrogen as the reference state, we were not able to find a 
stable structure for freestanding one-sided hydrogenation 
of graphane. This is no surprise, as the freestanding one- 
sided graphane is a highly polar structure. In order to 
find out whether the surface can balance the polarity and 
make the structure stable, we study the one-sided hydro- 
genation of graphene on the four different Si02 surfaces. 

First, the graphane hydrogen configurations on the 
four different Si02 surfaces were determined. One to six 
hydrogen atoms corresponding to 1/8-3/4 ML coverage 
of the graphene carbon atoms were placed in the super- 
cell on top of the carbon atoms. Then the carbon layer 
was deposited either on a top or hollow position with re- 
spect to the uppermost substrate atoms on a thinner slab 
corresponding to one unit cell of Si02 (5.3 A). The use 
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FIG. 1. The band structure along the T-K-K'-T direction, and the geometry of graphene on SiCb with different surface termina- 
tions. OHT: hydroxyl terminated, OT: oxygen terminated, SiT: silicon terminated, ROT: reconstructed oxygen terminated. In 
the band structures, the dashed line shows the band structure of freestanding graphene, and in the geometries the computational 
supercell. Atom colors: Black - carbon, blue - silicon, green - oxygen, yellow - hydrogen. 



of the thin slab allowed us to scan all different hydrogen 
initial configurations for both top and hollow positions, 
approximately 500 for each of the surfaces. The carbon 
layer and the hydrogens were relaxed. During the re- 
laxation, a large number of local energy minima were 
found and a tight convergence criterion proved to be ex- 
tremely important in order to locate the lowest-energy 
structures. Afterwards, the slab was increased to three 
Si02 unit cells (15.8 A) for the lowest-energy configu- 
rations on each surface and hydrogen coverage, and the 
three uppermost atom layers of the substrate were addi- 



tionally relaxed. For the calculation of the binding en- 
ergy, the reference state was chosen to be the optimized 
structure for graphene on the substrate, along with a cor- 
responding number of hydrogen molecules in a vacuum. 

Because the reconstruction on the ROT surface ex- 
tends deep to the substrate, we use the unhydrogenated 
optimized graphene- OT interface on the other side of the 
slab for the extensive configuration scan to decrease the 
computational effort. Due to the unsymmetry of the sur- 
faces, a dipole correction was used. Again, a large num- 
ber of initial hydrogen configurations was considered but 
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FIG. 2. (a) The binding energy (b) Distance of the carbons from the substrate (c) Band gap as a function of the hydrogen 
content for the lowest-energy configurations for the four surface terminations. In (b), the shaded regions show the distance 
range and the markers the average distance. 



during the relaxation of the carbon layer and the hydro- 
gen atoms, the highest-energy ones were discarded and 
relaxation was continued only for the low-energy config- 
urations. Later, the geometries were symmetrized and 
seven surface atom layers were relaxed. 

Fig. [2jt shows the binding energies of the lowest-energy 
configurations as a function of the hydrogen coverage. 
We see that a graphane layer with hydrogen only on one 
side has indeed stable structures on the Si02 surfaces, 
in contrast to the freestanding one-sided graphane. The 
surface termination, however, plays a profound role on 
the stability as there are no stable structures on the SiT 
and ROT surfaces. In general, the structures with an 
even number of hydrogen atoms per supercell are more 
stable than their odd- numbered counterparts, and we see 
an oscillatory behavior of the binding energy. This ap- 
plies also to the OHT surface if one takes into the account 
the hydrogen atom originating from the hydroxyl group 
of the substrate below the graphane layer. As there are 
multiple stable structures on the OHT surface, we note 
that even a small amount of hydrogen below the carbon 
layer stabilizes the structure and facilitates hydrogen ad- 
sorption above the carbon layer, also leading to higher 
hydrogen coverages. 

In general, unlike graphene, graphane is strongly 
bound to the substrate. This is illustrated in Fig. (2}d, 
where the average distance between the carbon and up- 
permost substrate atoms, 33 along with its variation, is 
shown as a function of the hydrogen content. On the 
reactive OT surface, the graphane layer is strongly cor- 
rugated and one stable structure (1/4 ML) is found. 
On the contrary, the unreactive ROT surface leads to 
long carbon-substrate distances and high binding ener- 
gies. The structures formed on the SiT surface are quite 
similar to the OT surface but not as strongly bound, as 
the average carbon-substrate distance is slightly longer, 
an exception being the 1/2 ML case where the distances 
are equal. On the OHT surface, the hydrogen atom below 
graphane facilitates the corrugation of the carbon layer 
and due to formation hydrogen-carbon bonds does not, 
in general, make the carbon-substrate distance longer. 33 



The band gap of graphane (Fig. [2k) does not increase 
monotonously when hydrogen coverage is increased but 
it shows a similar oscillatory behaviour as the energetic 
stability (Fig. [2^). The structures with large energy gaps 
are energetically more stable and have an even number of 
hydrogen atoms in the supercell (counting again one hy- 
drogen atom below the carbon layer on the OHT surface). 
Some of the hydrogen coverages show no energy gap at 
all. In most cases, this is due to relatively flat, almost 
nondispersive carbon bands at the Fermi level. It should 
be noted that DFT is known to underestimate band gaps 
and in the case of freestanding graphane, the use of the 
GW method leads to 1.5-2 times larger gaps j 19 i 2Q 

The three most stable structures on the OHT sur- 
face are, in the order of energetic stability, the 3/8, 5/8 
and 1/8 ML hydrogen coverages above the carbon layer. 
The band structures and hydrogen configurations of these 
structures are shown in Fig. [3ti-c. For the 1/8 ML cov- 
erage, the substrate hydrogen remains attached to the 
surface oxygen whereas for 3/8 and 5/8 ML coverages, 
the hydrogen atom binds from below to the carbon layer. 
The structure of these two states resembles that of free- 
standing graphane in the chair configuration as the hy- 
drogen atoms above the carbon plane bind only (mostly) 
to one sublattice for the 3/8 (5/8) ML structure, and 
most of the neighbours of these atoms bind either to the 
substrate or to a sub-carbon hydrogen. 

Also the OT surface forms a stable structure at 1/4 ML 
hydrogen coverage. In addition, the binding energy of 
graphane on the OT surface with 1/2 ML coverage is 
positive but close to zero, and thus the structure is pos- 
sibly a metastable state. Hydrogenated structures with 
1/4 ML coverage are lowest in energy also on SiT and 
ROT surfaces but their binding energies are significantly 
higher. 

On OT and SiT surfaces, the hydrogen atoms in 
the 1/4 ML configurations form rows along every other 
graphene zigzag line, the carbon atoms of the hy- 
drogenated lines alternatingly binding to a hydrogen 
atom above the carbon plane and to the substrate (see 
Fig. [3]ie). Interestingly, the corresponding band struc- 
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FIG. 3. Band structure and hydrogen configuration for the lowest-energy graphane structures. OHT: (a) 1/8 ML coverage (b) 
3/8 ML coverage (c) 5/8 ML coverage above the carbon layer. 1/4 ML coverage: (d) OT (e) SiT (f) ROT surface. Note that 
in the hydrogenated structures, the K and K' symmetry points are no longer equivalent. In (d) and (e), unhydrogenated zigzag 
carbon lines with a Peierls instability are formed (see text for details). Atom colors like in Fig. [1] 
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tures show a single Dirac cone between the K and K'- 
points of the first Brillouin zone, instead of the two cones 
seen for graphene at K and K\ A gap of magnitude 0.2 eV 
and 0.45 eV on the OT and SiT surfaces, respectively, ap- 
pears between the tips of the two cones , but the Fermi 
level is not within this gap, reducing the indirect gap to 
few meV. The lowest-energy hydrogenated structure on 
the ROT surface, on the contrary to the OT and SiT 
counterparts, is a semiconductor with a 2.5 eV band gap 
and does not show any significant directional variation in 
the band structure. The hydrogen atoms adsorb only to 
one sublattice but they do not form lines. 

We gain more insight into the merging of the two Dirac 
cones of 1/4 ML OT and SiT structures by plotting the 
electronic bands both along the direction of the hydro- 
genated chains, and perpendicular to it. In the trans- 
verse direction, the bands are almost flat and nondis- 
persive whereas along the direction of the zigzag chains, 
the Dirac cone is present. Effectively, the systems are 
thus one-dimensional conductors, although some cou- 
pling over the hydrogen lines through the substrate re- 
mains. Very recently,— a Peierls instability was found 
in similar unhydrogenated zigzag chains in freestanding 
graphane. The instability is seen also in our structure, 
with the bond length alternation along the zigzag chain 
being of somewhat smaller amplitude than in the free- 
standing structure (Ad = 0.008-0.012 A compared to 
0.015 A 3 - 5 ). This is expected, as the bond lengths of 
graphane on the substrate are generally compressed in 
comparison to the freestanding structure. 

The two graphene Dirac cones have been suggested 
to be useful for valley tronics^ an analog of spintronics 
using the valley isospin arising from the inequivalence 
of the two Dirac cones. Topological line defects have 
been suggested to give rise to valley filtering. 37 The 
valley degree of freedom is also thought to be robust 
against disorder^ As there is only a single valley in 
the hydrogenated structure, our results suggest that the 
valley degree of freedom could be destroyed by ordered 
adsorption. 



V. THERMODYNAMIC AL STABILITY 

In contrast to the experiments, our calculations are 
performed in vacuum at K. The experimental condi- 
tions, such as the hydrogen content of the environment, 
naturally affect the outcome of the hydrogenation. Using 
an ab initio thermodynamics approach ^ 38 ! 39 we may ad- 
dress the relative stability of the different hydrogenated 
surfaces as a function of the chemical potential of the 
hydrogen species. The Gibbs energy of the hydrogen ad- 
sorption is calculated as 
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FIG. 4. (a) The Gibbs energy of adsorption on the four sur- 
faces as a function of the chemical potential of molecular hy- 
drogen (lower x-axis) and atomic hydrogen (upper x-axis). 
For positive AG a ds, the unhydrogenated surface is favored. 
The preferred structures with the lowest AG a ds are marked 
with solid lines and the corresponding coverage as monolay- 
ers is shown, (b) (Color online) The chemical potential of 
molecular hydrogen H2 (blue/black) and atomic hydrogen 
(green/gray) in the gas phase as a function of temperature 
at three different partial pressures of the hydrogen species. 



tion of adsorbed hydrogen (in the units of hydrogen 
atoms per graphene carbon) and Ph 2 /h — ^°h 2 /h 

(T) + 

fceTm Pli pi H is the chemical potential of hydrogen 
molecules or atomic hydrogen, where P is the partial 
pressure of the hydrogen species. The standard chemical 
potential M# 2 (T) can be obtained from thermodynami- 
cal reference data— and is shown in Fig.0]3. We can also 
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compare the hydrogenation through hydrogen molecules 
to atomic hydrogen by defining the binding energy and 
AG a ds with respect to atomic hydrogen in gas phase. 

AGads as a function of hydrogen chemical potential for 
each of the four surfaces is shown in Fig. HJ considering 
both molecular (H2) and atomic (H) hydrogen. The 
structure with lowest AG a ds is the most favorable. At 
high chemical potentials, the structures with a higher 
hydrogen content become favored, even if their DFT 
binding energies are positive and the structures are thus 
unstable in vacuum at K. Due to the logarithmic 
dependence of the chemical potential on the pressure, 
reaching positive values for the chemical requires high 
temperatures and pressures. Thus, the more densely 
hydrogenated structures are hard to reach in experi- 
ments using molecular hydrogen. The stability diagrams 
presented in Fig. |4] are, of course, idealistic and do not 
take into account thermal kinetics but they could serve 
as a rough guideline to the experimental conditions 
required to achieve different hydrogen coverages. 



VI. CONCLUSIONS 

In conclusion, we have studied graphene and its hy- 
drogenated version, graphane, on the a-quartz Si02 
substrate, considering four different terminations of the 
(0001) surface. In our calculations, graphene does not 
form bonds with the substrate atoms. Instead, the linear 



dispersion in the vicinity of the K and K' points of the 
first Brillouin zone is retained. 

We find that the surface termination has a profound 
effect on the hydrogen adsorption. On the OHT and OT 
surfaces, stable coverages are found, with 3/8, 5/8 and 
1/8 ML, and 1/4 ML hydrogen content, respectively. The 
presence of hydrogen below the carbon plane facilitates 
the adsorption, stabilizing the hydrogenated structures 
and allowing higher surface coverages. The higher cover- 
ages (3/8 and 5/8 ML) on the OHT surface are semicon- 
ductors with band gaps 3.5-3.9 eV, whereas on the OT 
surface at 1/4 ML coverage we observe a peculiar band 
structure with a single Dirac cone between the K and K' 
points. 

Our results indicate that in order to hydrogenate 
graphene deposited on Si02, it is beneficial to oxy- 
genate the surface and passivate it with hydrogen be- 
fore graphene deposition. As Si02 is amorphous, uneven 
hydrogen adsorption may be explained through surface 
domains of different termination. In addition, a small 
amount of hydrogen under the graphene layer is not likely 
to change the conducting properties of graphene. This 
means that partial dehydrogenation 7 of graphane on the 
hydrogenated surface is still an option and could be useful 
in nanoelectronics, for instance for forming nanoribbons 
in an insulating graphane matrix i 41 i 42 
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